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A theoretical model for the phonon-mediated off-resonant coupling between a quantum dot and 
a cavity, under resonant excitation of the quantum dot, is presented. We show that the coupling 
is caused by electron-phonon interaction in the quantum dot and is enhanced by the cavity. We 
analyze recently observed resonant quantum dot spectroscopic data by our theoretical model. 



One of the interesting recent developments in cav- 
ity quantum electrodynamics (CQED) experiments with 
quantum dots (QDs) coupled to semiconductor micro- 
cavities is the observation of off-resonant dot-cavity cou- 
pling. This unusual phenomenon is observed both in 
photoluminescence studies, where an above-band laser 
generates carriers that incoherently relax to rccombinc 
at the QD frequency [1-3] and also under resonant ex- 
citation of the QD or the cavity [4, 5]. The coupling 
observed via photoluminescence is attributed to sev- 
eral phenomena including pure QD dephasing [6], the 
electron-phonon interaction [7, 8], multi-exciton com- 
plexes [9] and charges in the vicinity of the QD [10]. To 
isolate the role of phonons in off-resonant QD-cavity cou- 
pling, studies employing resonant excitation of the QD 
are preferable as they avoid possible complications aris- 
ing from multi-excitonic complexes and nearby charges 
generated in above band pumping. Resonant excitation 
of a QD coupled to an off-resonant cavity has been re- 
cently used to perform QD spectroscopy, enabling the 
observation of power broadening of the QD line-width 
and saturation of the cavity emission [11-13]. However, 
there is presently no theoretical model accounting for off- 
resonant dot-cavity coupling for the case of resonant QD 
excitation. Although off-resonant coupling can be mod- 
eled by introducing a phenomenological incoherent cav- 
ity pumping mechanism [14], such a treatment masks the 
actual physical phenomenon responsible for the coupling. 
Without an explicit coherent driving term in the system 
Hamiltonian, the resonant QD spectroscopy results can- 
not be explained. 

In this Letter, we theoretically model the phonon- 
mediated interaction between a cavity mode and an off- 
resonant QD under resonant excitation of the QD. We 
first model the coupling via pure QD dephasing. Then, 
we propose a new model where the phonon-mediated cou- 
pling is enhanced by the presence of the cavity. We com- 
pare these two models and find that they provide quali- 
tatively similar signatures in terms of experimental reso- 
nant QD spectroscopic studies, such as power broadening 
and saturation of a resonantly driven QD [11, 12]. How- 
ever, in the newly proposed model, the coupling is main- 
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tained at very large QD-cavity detunings (~ 3 meV), as 
observed in the recent experiments [11]. We also observe 
the signature of the inherent asymmetry between phonon 
emission and absorption rate in the simulated QD spec- 
troscopy results, depending on whether the QD is red or 
blue detuned from the cavity. 

The dynamics of a coherently driven QD (with a 
ground state \g) and an excited state |e)) coupled to an 
off-resonant cavity mode is governed by the Hamiltonian 
H (in a frame rotating at the driving laser frequency) 

H = Aoj^a + AtUaO-^a + ig(a)a - aa f ) + 0(cr + cr 1 *) (1) 

where, a and a are the annihilation and lowering op- 
erators for the cavity mode and the QD, respectively; 
Alo c = uj c — loi and Aw a = u a — uji are the cavity 
and dot detunings from the driving laser, respectively; 
A = uj a — uj c is the QD-cavity detuning; is the Rabi 
frequency of the driving laser and is proportional to the 
laser field amplitude and g is the coherent interaction 
strength between the QD and the cavity. 

In this coupled system, there are two independent 
mechanisms for energy dissipation: cavity decay and QD 
dipole decay. The system losses can be modeled by the 
Liouvillian, and the Master equation describing the dy- 
namics of the lossy system is given by 

^ = -i[H,p]+2KC[a] + 2 1 L[a] (2) 

where p is the density matrix of the coupled QD-cavity 
system, 27 and 2k arc the QD spontaneous emission rate 
and the cavity population decay rate, respectively. We 
neglect any non-radiative decay of the QD exciton. C[D] 
is the Lindblad operator corresponding to a collapse op- 
erator D and is given by: 

C[D] = Dprf - \^Dp - X - P D^D (3) 

In addition, phonons in the solid state system destroy 
the coherence of the exciton. This is generally modeled 
by adding an additional incoherent decay term 2^d£. [o~^ <r\ 
to the Master equation, where 27^ is the pure dephasing 
rate of the QD. This term destroys the polarization of 
the QD without affecting the population of the QD. The 
dissipation of the QD polarization and population (a z = 



2 



[crt, c]) is given by the following mean- field equations: 



virtual(3) 



dt 
dt 



= -(7 + 7d)(cr) 
= -2 7 (l + (£r,» 



(4) 
(5) 



Hence, the linewidth of the QD, at the zero excitation 
power limit, is given by 2(7+7^). However, in this model, 
the effect of phonons is embedded in the phenomcnolog- 
ical pure dephasing rate 7^, which affects only the QD, 
and does not include any cavity effects. 

Wc now propose a different model for off-resonant dot- 
cavity coupling, where the phonon-mcdiatcd coupling 
strength is affected by both the cavity and the QD. 
The effect of phonons can be modeled by replacing the 
pure dephasing term with two additional incoherent de- 
cay terms in the Master equation. For blue-detuned QD 
(Fig. f a) we have: 



2n£l 



+ 27 r n£(cr t a) + 2^ r (n + 1)£(W) 



(6) 



where, 27 r is an effective decay rate of the QD exci- 
ton states via the emission of a phonon and a photon 
at the off-resonant cavity frequency, n is the average 
number of phonons at the dot-cavity detuning frequency 
(A) present in the system at thermal equilibrium with 
the reservoir at a temperature T, and is given by 



n(A,T) 



f 



e hA/k B T _ 1 



(7) 



The analysis for a QD red detuned from the cavity (Fig. 
1 b) can be carried out in a similar manner by replacing 
the final two terms of Eqn.6 with 2^ r nL(aa^) and 2 7r .(n+ 
l)£(<jta). 

The decay term a denotes the annihilation of a cav- 
ity photon and excitation of the QD, while the term era' 
denotes the creation of a cavity photon and collapse of 
the QD to its ground state, accompanied by the creation 
(or annihilation) of phonons to compensate for the QD- 
cavity frequency difference. Only the second process is 
important for observing cavity emission under resonant 
excitation of the dot. We also note that the observation 
of QD emission under resonant excitation of the cavity 
[4] can be modeled in the same way by changing the co- 
herent driving term from fl(er + a^) to Q(a + a'). In this 
situation, the collapse operator er^a is important. Simi- 
lar decay channels have been proposed to model cavity 
assisted atomic decay [15]. A detailed derivation of these 
incoherent terms can be found in the Supplementary Ma- 
terials. We call this process a cavity-enhanced phonon 
process. 

We first consider the case where the QD is blue detuned 
from the cavity (Fig. 1 a). The qualitative nature of the 
dissipation of the QD polarization and population can be 
determined by the mean-field equations (assuming n = 0, 
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FIG. 1: (a),(b) Level diagram of the coupled QD/cavity sys- 
tem for blue (a) and red (b) detuned QD relative to the cav- 
ity resonance. A laser drives the quantum dot (transition 
between ground state (1) and excited state (2)) resonantly. 
The excited state (2) can decay via two paths: the first is 
by direct decay back to the ground state (1) via the sponta- 
neous emission of the QD; the second is by indirect decay via 
the emission (Fig. la) or absorption (Fig. lb) of a phonon 
(transition (2) to (3)) and subsequent emission of a photon at 
the cavity frequency (transition (3) to (1)). (c), (d) Emission 
S(ui) as a function of frequency uj for a QD resonantly driven 
by a laser (Acj a = 0). The plot (c) corresponds to a QD that 
is blue detuned, and (d) corresponds to a QD that is red de- 
tuned from the cavity. Three different cases are considered: 
first, without any dephasing; second with pure dephasing and 
third with a cavity enhanced phonon process (our newly in- 
troduced model). We observe the Mollow triplet at the QD 
frequency in all three cases. The inset of (c) and (d) show an 
enlarged view of the emission at the cavity frequency. We ob- 
serve no off-resonant cavity emission without pure dephasing, 
or our newly introduced cavity enhanced phonon process. For 
the simulation we assume g/2n = n/2n = 20 GHz; 7/2-7T = 1 
GHz; Q/2tt = 6 GHz; QD-cavity detuning A/2tt = ±6 GHz. 



i.e., at the zero temperature limit): 
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We notice that, unlike the pure dephasing case, both 
the QD population and polarization are affected by 
the cavity enhanced phonon process. The linewidth of 
the QD at the zero excitation power limit is given by 
27 + 277.(1 + (a^a)), which is qualitatively different from 
what we expect from a pure QD dephasing model (as 
the presence of cavity photons affects the QD linewidth). 
We emphasize that for both models, the QD lincwidths 
obtained from the mean-field equations neglect the modi- 
fication resulting from cavity QED effects in the presence 
of the cavity and hence the predicted linewidths are only 
recovered at large QD-cavity detunings. 

The resonance fluorescence of the system is given by 
the power spectral density (PSD) of the coupled system. 
As we collect the fluorescence primarily from the cavity, 
the PSD is calculated as the Fourier transform of the 
cavity field auto-corrclation: 

/oo 
(a t (r)a(0))e- MT dr (10) 
-00 

To determine the two time correlation functions and, sub- 
sequently, the PSD, we use the quantum regression the- 
orem [16]. 

We simulate the coupled system using numerical inte- 
gration routines provided in the quantum optics toolbox 
[17] with realistic system parameters k/2it = g/2n = 20 
GHz and 7/27r = 1 GHz. The off-resonant coupling is ob- 
served both for a strongly (g > k) and a weakly coupled 
QD to the cavity (g < k). Figure 1 c shows the numeri- 
cally calculated resonance fluorescence spectra obtained 
from the cavity under resonant excitation of a blue de- 
tuned QD for three different cases: no dephasing, pure 
QD dephasing and phonon-induced exciton decay at a 
bath temperature of 4 K. We first note that no emission 
is observed at the cavity frequency in the absence of pure 
dephasing or cavity enhanced phonon process (see inset 
of Fig. 1 c). Though the pure dephasing (■yd/2ir = 1 
GHz) and phonon-induced decay (j r /2TT = 1 GHz) cases 
both show off-resonant cavity emission, the latter shows 
enhanced cavity emission. In all three cases, we also ob- 
serve the Mollow triplet at the QD frequency, as expected 
from QD resonance fluorescence [18]. The Mollow side- 
band closer to the cavity is enhanced compared to the 
other side-band causing an asymmetric triplet. Fig. 1 
d shows the similar spectra for a red-detuned QD. Our 
model predicts that the emission at cavity frequency is 
considerably smaller for a red detuned QD compared to 
a blue-detuned one as expected at finite temperature (as 
the process for a red-detuned dot relies on phonon ab- 
sorption, while for a blue-detuned dot relies on phonon 
emission). 

We then perform simulations to observe how the emis- 
sion collected at the cavity frequency depends on the 
dephasing and phonon-induced decay rates, as well as 
on the field of the driving laser resonant with the QD 
(Auj a = 0). We observe that the cavity emission I in- 
creases almost linearly with the pure QD dephasing rate 
7<j [Fig. 2 a] , but exhibits a nonlinear dependence on the 



rate "f r when the coupling is enhanced by the presence of 
the cavity. Fig. 2 b shows the cavity fluorescence as a 
function of the laser Rabi frequency fi. For both models, 
the cavity fluorescence / follows a saturation curve 



where, P oc fl 2 and I sa t is the saturated cavity emission 
intensity. As noted previously in the article, the cavity 
emission is higher when the process is enhanced by the 
presence of the cavity. Furthermore, we investigate the 
dependence of the saturation emission intensity I sat as a 
function of the QD-cavity detuning A (Fig. 2 c). I sa t 
falls of as 1/ A 2 with the detuning A when the dot-cavity 
coupling is modeled as a pure dephasing process. How- 
ever, when the coupling is modeled as a cavity enhanced 
phonon process, the saturation intensity exhibits a di- 
minished dependence on detuning A (estimated to be 
~ 1/A 25 ). Hence, one may observe off-resonant cou- 
pling for larger detunings when the phonon process is 
enhanced by the cavity. Fig. 2 d shows log(I) as a func- 
tion of the cavity decay rate n, for a fixed detuning of 
A/2-7T = 200 GHz. For both models, the emission falls 
off as 1/k 2 , signifying that the off-resonant coupling does 
not depend on the overlap between the QD and the cavity 
spectra. 

We now measure the QD line-width Auj monitoring 
the cavity emission, while scanning the laser wavelength 
across the QD resonance, similar to the experiments 
in Ref. [11]. We observe that at very low excitation 
power f2/27r = 1 GHz and large QD-cavity detuning 
A/27T = 12k, the lincwidths of the QD are very close 
to the theoretical linewidth in the absence of the cavity 
(shown by the solid black line in Fig. 3 a). At a constant 
QD-cavity detuning and laser excitation power, the QD 
line- width increases with increasing 7^ and 7 r (Fig. 3 a). 
We observe broadening of the line-width with increas- 
ing laser power (Figs. 3 b). The power broadened QD 

linewidth Auj is fit with the model Auj = AujqVI + P, 
where Aujq is the intrinsic line-width of the QD and P 
is obtained from the fit to the saturation of the cavity 
emission [11]. The theoretical model does not reproduce 
the additional power-independent broadening of the QD 
[11] and we believe that this extra broadening may result 
from QD spectral diffusion [19]. We analyze the intrin- 
sic QD linewidth Aujq (without power broadening, i.e., 
obtained from plots in Figs. 3 b at SI = limit) as a 
function of the dot-cavity detuning A for two different 
models (Figure 3 c). We observe that at large detuning, 
Aujq approaches the unperturbed QD linewidth 2(7 + 7^) 
and 2(7 + 7,,), respectively. We fit empirical models of 
A~ a to the intrinsic linewidths for the two models and 
find that with pure QD dephasing, Awo falls off more 
slowly (a ~ 0.4) compared to the cavity enhanced cou- 
pling (a ~ 0.7). The weak dependence of the intrinsic 
QD linewidths on the dot-cavity detuning shows that the 
off-resonant cavity does not perturb the QD significantly. 
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FIG. 2: Resonance fluorescence / collected from the cavity for 
a blue-detuned QD (A > 0). (a) log(I) as a function of the 
rates jd and y r (for two models, respectively). A/k = 10. (b) 
Normalized cavity fluorescence as a function of the Rabi fre- 
quency of the laser for two models. Saturation of the cav- 
ity emission is observed. The fluorescence values are not in 
scale, (c) Dependence of the saturated cavity emission I sa t 
on the dot-cavity detuning A. (d) log(I) as a function of 
the cavity linewidth k. (Parameters used for all the simu- 
lations are g/2-ir = k/2tv = 20 GHz; y/2-ir = 1 GHz. For 
Figs. 2 a, c and d, Q,/2tt — 4 GHz and for Figs. 2 b,c and d, 
7d/27r = 7 r /27r = 1 GHz.) 



However, the results are dramatically different for a 
red-detuned QD. This is a result of the fact that the rate 
of the incoherent process involving the operator aa^ is 
different depending on whether the QD is red or blue de- 
tuned from the cavity (because at any temperature, the 
rates of absorption and emission of phonons arc differ- 
ent). This asymmetry is manifested in both the emission 
collected at the cavity resonance and in the QD spec- 
troscopy result. Fig. 4 a shows the difference in the QD 
linewidths as a function of the bath temperature T for 
different driving laser Rabi frequencies f2, when the QD 
is red or blue detuned from the cavity by same abso- 
lute value. We observe that the difference in linewidth 
is higher when the QD is weakly driven and hence is not 
power broadened. Full quantum optical simulations show 
that the linewidth depends on the bath temperature, but 
reaches the maximum value of 2j r at a higher tempera- 
ture. Fig. 4 b shows the ratio of the cavity emission as 
a function of the bath temperature for different driving 
laser Rabi frequencies tt. The difference in cavity inten- 
sity is maximum at lower bath temperature and is almost 
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FIG. 3: Linewidth of the blue detuned QD relative to the cav- 
ity (A > 0) as measured by monitoring the off-resonant cavity 
emission (similar to the experiments in Ref. [11]). Results ob- 
tained from two models are presented, (a) QD linewidth as 
a function of the rates 7^ and ~y r for pure dephasing and the 
cavity enhanced phonon process, respectively. In both cases, 
Q,/2n — 1 GHz and A = 12k. The solid black line shows 
the theoretical estimates of the QD linewidth when the laser 
excitation power is very low and the QD is not significantly 
perturbed by the cavity, (b) Power broadened linewidth vs 
laser Rabi frequency for QD-cavity detuning A/2n = 6 
GHz (for the case of pure QD dephasing and the cavity en- 
hanced phonon process), (c) Dependence of the intrinsic QD 
linewidth Auo on dot-cavity detuning A. Aujo/2it approaches 
2(7 + 7d)/27r or 2(7 + -y r )/2n ((both chosen to be 4 GHz) 
with large A. (Parameters used for all the simulations are: 
k/2tv = g/2n = 20 GHz; 7/2tt = 1 GHz.) 



zero at higher temperature. 

In conclusion, we have proposed a theoretical model for 
the off-resonant dot-cavity coupling under resonant exci- 
tation of the QD. Introduction of an incoherent decay 
channel (referred to as the cavity enhanced phonon pro- 
cess, which is different from pure dephasing) shows that 
the phonon-mediatcd dot-cavity coupling is enhanced by 
the presence of the cavity. By comparing the power- 
broadening and saturation of the QD between the cav- 
ity enhanced phonon process and pure QD dephasing, 
we found that the dot-cavity coupling is enhanced in the 
former case and is observed for a larger QD-cavity detun- 
ing. Our model can also be used to explain asymmetry in 
the spectroscopy results for a QD blue and red detuned 
from the cavity. We believe that such an off-resonant 
dot-cavity coupling can be used as an efficient read-out 
canncl for resonant QD spectroscopy and for QD-spin 
manipulation. 

The authors acknowledge financial support provided 
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FIG. 4: (a) The difference in QD linewidths measured via col- 
lected emission through off-resonant cavity, for a blue (Au+) 
and a red (Auj- ) detuned QD for different values of the excita- 
tion laser Rabi frequency Q as a function of the bath temper- 
ature T. (b) Ratio of the cavity intensity for an off-resonant 
QD blue (7+) and red detuned from the cavity, as a func- 
tion of the bath temperature T. For all the simulations, the 
absolute value of QD-cavity detuning is kept at 10/t. 



transition and g 3 23 is the interaction strength between the 
QD exciton and an j th mode phonon. We note that this 
interaction Hamiltonian is valid only for the level struc- 
ture as in Fig. 1 a, where the cavity is at lower energy 
than the QD. For the situation in Fig. 1 b (where the 
cavity is of higher energy compared to the QD), the in- 
teraction Hamiltonian Hi is given by 

F 7 = ^(a|3)(l| + a+|l)(3|)+^ fl i 3 (6 J -|3)(2| + 6]|2)(3|) 

(15) 

In the following derivation, we will use the situation 
shown in Fig. 1 a. 

If we define the QD resonance frequency as uj ai then 
u a = ^2 — <^x; and the cavity frequency is given by u c = 
lo 3 — u)i . Then the QD-cavity detuning is given by A = 
LO2 — ^3. Defining a\j — we can write 

ois = -i[vi3,H + Hi] 

= -iuj c <r 13 - ig v a(an - a 33 ) - i £ g 3 23 b\(Ji2 

j 

Similarly, 

023 = -i[(T23,Ho + Hi] 

= iA<r 23 - ig v aa 21 - i £ fli^C ^ - 033) 
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I. SUPPLEMENTARY MATERIALS 

A. Derivation of Decay Terms 

We use the level diagram as shown in Fig. 1 a to model 
the effect of phonons explicitly. The Hamiltonian of the 
system is given by 



H = H + Hi 



(12) 



where, 



Ho = wi|l)(l| +«a|2)<2| +w 3 |3)(3| + woW J^i/jfcJ&j 

(13) 



and 

F 7 = ft,(o|3)<l| + o+|l)(3|) 



£s^]|3><2|+^|2)(3|) 



where, is the population operator for i level 



(14) 
a is 



the annihilation operator for the cavity mode; bj is the 
annihilation operator for a phonon in the j mode. u>i , 
u> and Vj are the frequencies of the i th energy level, cavity 
resonance and a phonon in the j th mode. g v signifies the 
interaction strength between the cavity and the virtual 



Separating the slow and the fast components of the op- 
erators, we can write 



013 
023 
012 

a 



ai 3 e 
023 e 



-ioj c t 

-ibJjt 



<7i2e 



p e -io>jt 



i(ui c —ujj)t 



(16) 
(17) 
(18) 
(19) 
(20) 



Hence the equations governing the dynamics of the sys- 
tem can be written as 

013 = -ig v a(an - <t 33 ) - i V] g^M&i-i (21) 



3 



and 

O23 = l(A — Wj)cT23 — i5t)2cr2i- 



^523^(^22-033) (22) 



As level 3 is a virtual level, it is never populated. Hence 
by adiabatic elimination, using 013 = 023 = 0, we obtain 



023 



g v ad 2 i + Ej .923^(^22 - 0-33) 
A — ujj 



and 



012 = 



-g v a(a n 



033 J 



Y, 3 9i 3 b) 



(23) 



(24) 



6 



Using these values, we can find the interaction Hamilto- 
nian. The first term <7„( acr 3i + o 013) denotes the coher- 
ent dynamics. The second term, which signifies the effect 
of phonons, can be written as (using the Eqns. 23 and 
24) 



and 



H ph = ^2gi 3 (b]a-32 + bj<T23) 
3 

3 



(25) 
(26) 

(27) 



E 



(fa-twXbm + bjbj) (28) 



The second term involves two-phonon processes which 
are less likely. If we neglect them, we can model the 
effect of phonons as follows: 



H, 



(Ma t <Ti2 + 6j5<72i) 



We can write this Hamiltonian as 

H p h = (3*1712? f + OCT 2 lf ) 

where, the operator T can be written as 



and 



E 



9239v 



bje- 



(29) 



(30) 



(31) 



(32) 



To obtain the familiar Lindblad term, we take the partial 
trace of the correlation between the reservoir operators 
over the reservoir variables. The correlation is given by 



<r t (f)r(t))fl = E 



j 



A — uj.i 



(33) 



As the operators 6j are bosonic and the system is in ther- 
mal equilibrium with a bath at temperature T, using the 
relation 



(b]bj) R = n(u)j,T) 



1 



(34) 



- 1 



we find 



3233f 



A — w, 



-<«iC*-t') fi ( a ,. j T) (35) 



<r(Ort(t)> fl = X; 



9239" 



-*"jlt-t)(n( Ujt T) + i) 
(36) 

Prom the correlation, we find that the phonons with fre- 
quency A (corresponding to the difference between lev- 
els |2) and |3), i.e., QD-cavity mode detuning), have the 
maximum contribution in the interaction Hamiltonian. 
In the Born-Markov approximation, we can model the 
electron-phonon interaction (for Fig. 1 a) as an inco- 
herent decay process by adding two extra terms to the 
Master equation: 2j r nC(a^a) and 2-f r (h + l)£(cr(v), j r 
being the effective decay rate of the excited QD state and 
is given by 



It 



32 3 3» 



A-, 



(37) 



For the situation shown in Fig. 1 b, the decay terms are 
given by 2j r hC(aa ji ) and 2'y r (n+l)C(a' a). We note that 
the different rates in both cases are due to an inherent 
asymmetry between the absorption and emission rates of 
the phonons. 



B. Derivation of the Mean Field Equations 

To find the mean field equations for an operator A from 
the Master equation, we used the following relation: 



d(A) 
dt 



jTr[Ap] = Tr 



dt 



(38) 



For the cavity enhanced phonon process, the mean field 
equations for a non-zero (n) is given by (when the QD is 
blue detuned from the cavity) 

d(a) 



dt 
dt 



= -j(<j)-j r (l + (^a))(a) (39) 

7 r n(l + 2(a t a»(cr} (40) 

= -2 7 (1 + (<j z }) (41) 

-2 7r (l + n)(l + (at a ))(l + (a z }) (42) 



When the QD is red detuned from the cavity, the mean 
field equations are: 



d{°) 
dt 

d(<r z ) 
dt 



-7(d) - 7 ? ,(a t a)(er) 
7 r n(l + 2(a t a))(er) 
-2 7 (1 + (a z )) 
-27 r n(l + (a f a})(l + (a z )) 



(43) 
(44) 
(45) 
(46) 



We note that while deriving these mean-field equa- 
tions, we assume that the cavity and QD operators are 
uncorrclatcd and write 



{a) aa) = («ra)(cr) 



(47) 
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Under weak driving, the approximation holds very well. make any assumptions. 
However, in full quantum optical simulations, we do not 
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